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We analyze van der Waals interactions between two rigid polymers with sequence-specific, 
anisotropic polarizabilities along the polymer backbones, so that the dipole moments fluctuate 
parallel to the polymer backbones. Assuming that each polymer has a quenched-in polarizability se¬ 
quence which reflects, for example, the polynucleotide sequence of a double-stranded DNA molecule, 
we study the van der Waals interaction energy between a pair of such polymers with rod-like struc¬ 
ture for the cases where their respective polarizability sequences are (i) distinct and (ii) identical, 
with both zero and non-zero correlation length of the polarizability correlator along the polymer 
backbones in the latter case. For identical polymers, we find a novel scaling behavior of the van 
der Waals interaction energy for small inter-polymer separation r, in contradistinction to the r~* 
scaling behavior of distinct polymers, with furthermore a pronounced angular dependence favoring 
attraction between sufBciently aligned identical polymers. Such behavior can assist the molecular 
recognition between polymers. 

PACS numbers: 82.70.Dd, 83.80.Hj, 82.45.Gj, 52.25.Kn 


I. INTRODUCTION 

Van der Waals (vdW) forces are ubiquitous in Nature, 
being especially important in the nanoscale world [1]. 
They arise owing to the presence of fluctuating perma¬ 
nent and/or induced dipoles on atoms, molecules and 
molecular aggregates. Dipolar fluctuations enter the cal¬ 
culation of the vdW interaction (free) energies through 
the corresponding frequency-dependent and, in general, 
anisotropic dielectric response functions PHD ■ The latter 
case is particularly germane for long chain-like molecules, 
where the dielectric response functions along the molec¬ 
ular chain axis and perpendicular to it can differ signif¬ 
icantly. A good example would be the strong optical 
anisotropy of DNA [5] or long carbon nanotubes [6] with 
significantly different dielectric response functions in the 
axial and the radial directions. The role of dielectric 
anisotropy in the context of polymers was recently rec¬ 
ognized specifically in the context of liquid crystalline or¬ 
dering in the case of anisotropic polarizable polymers [7] , 
leading to important consequences. In general, vdW in¬ 
teractions between anisotropic media lead not only to 
long-range interactions but also to long-range torques as 
has been recognized some time ago [iH]. The details of 
this anisotropy effect are in fact many faceted and are 
discernible in the non-retarded as well as in the retarded 
limit of vdW interactions 

Anisotropy is not the only important peculiarity of 
polymer polarizability. It is now recognized, based 
on extensive ab initio calculations and UV-vis molar 
absorbance measurements, that the specificity of the 
oligonucleotide sequence of DNA molecules strongly in¬ 
fluences the dielectric response of the molecule and its 
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polarizability properties m- The question thus arises 
whether the detailed sequence of monomer polarizabili¬ 
ties in heteropolymers would influence in some essential 
capacity the interaction between two identical polymers, 
enabling in this way the interacting biological systems 
to recognize at a finite separation the polarizability pat¬ 
tern, leading to molecular recognition. Similar types of 
molecular recognition problems have been studied in the 
context of electrostatic interactions between randomly, 
irreversibly charged objects [TTHT^ and randomly pat¬ 
terned surfaces interacting through a generic finite range 
interaction potential, as a prototype model of protein- 
protein interactions [Blin]. It seems, however, that at 
least for polynucleotides, the sequence specificity is much 
more pronounced in terms of polarizability than charge 
distribution and thus there is some interest to investi¬ 
gate the sequence-specific polarizability effects on vdW 
interactions. 

For our purpose we consider a pair of rigid, charge 
neutral and polarizable polymers (or, more precisely, het¬ 
eropolymers), which experience induced dipole-induced 
dipole interactions between monomers characterized by 
different values of polarizability. In general, the polariz¬ 
ability is a tensor, with different radial and axial compo¬ 
nents (the axial direction being parallel to the backbone 
of the polymer). For long and thin polymers, the polar¬ 
ization of the polymer is mostly concentrated along its 
backbone, and one can then neglect the radial contribu¬ 
tion to the polarizability, which is the approximation that 
we will consider in this paper. The local polarizability of 
a heteropolymer depends on the identity of the particular 
monomer (e.g., base pairs along the double-stranded (ds) 
DNA molecule), and we can thus take it to be set irre¬ 
versibly, i.e., quenched. For simplicity, we consider only 
the case of non-retarded interactions since the molecu¬ 
lar recognition effects are expected to operate mainly on 
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short length scales. Furthermore, we consider the regime 
of high temperatures and focus on the contribution of 
static polarizabilities, which are nevertheless anisotropic 
so that the molecular dipoles can not be considered as 
freely rotating [5T]. Quite generally, in aqueous environ¬ 
ment the static term contributes more then 50 % to the 
total vdW interaction potential [1]. Finally, we delimit 
our analysis to the case of rod-like polymers, which we 
will study using the so-called shish-kebab model [52] ; this 
model should also be appropriate for dsDNA molecules 
at small separations. 

We shall derive a field theory that accounts for ther¬ 
mal fluctuations of induced dipoles along the polymer 
backbone, which we shall then apply to study the effec¬ 
tive interaction energy between a pair of rod-like poly¬ 
mers with distinct (i.e., completely uncorrelated) polar¬ 
izability sequences or a pair of polymers with identical 
(i.e., maximally correlated) polarizability sequences. By 
comparing the vdW interaction energies of a pair of dis¬ 
tinct polymers and a pair of identical polymers, we find 
the surprising result that the interaction energies dis¬ 
play fundamentally different behaviors when the inter¬ 
polymer separation r is smaller than the length of each 
polymer. T he se quence-averaged (in the sense to be de¬ 
fined in Sec. IIB) interaction energy for a pair of distinct 
polymers shows a dependence and is invariant with 
respect to the inversion of either polymer about its cen¬ 
ter of mass, whereas the sequence-averaged interaction 
energy for two identical polymers has a scaling form 
and is non-degenerate with respect to polymer inversion. 
The latter effect leads to favoring the selection of pairs of 
identical over distinct rod-like polymers and, moreover, 
amongst those pairs of identical polymers selected, there 
is a higher probability that these polymers are (i) parallel 
or (ii) anti-parallel with the line joining their centers of 
mass parallel to the tangent vector of each. Such features 
can assist in the recognition of molecular sequences, such 
as in DNA molecules, interacting at short distances. 

Here is a brief outline of the paper. In Sec. II, we 
describe our model of rigid, polarizable and charge neu¬ 
tral polymers interacting via induced dipoles. We make 
an interesting observation about the character of interac¬ 
tions between induced dipoles along charge neutral semi- 
flexible polymers, and present a discussion of sequence- 
specific polarizabilities, describing two models for the 
correlation between such polarizabilities. In Sec. Ill, we 
derive a field theory for the interactions of polymers via 
induced dipoles fluctuations. In Sec. IV, we present our 
results for the effective interaction energy between a pair 
of rigid, rod-like polymers. In Sec. V, we give a summary 
and discussion of the key results. 


II. MODEL DESCRIPTION OF RIGID 
POLARIZABLE POLYMERS 


dsDNA molecule) on the f-th polymer by the symbol 
(to). We assume that the polymer has M monomers 
(so TO = 1,..., M) each of equal length a, and there are 
a total of N polymers (so i = 1,... ,N). Each segment 
comes equipped with a thermally fluctuating, instanta¬ 
neous dipole moment represented by a vector (to) and 
a polarizability We assume that the dis¬ 

tribution along the polymer backbone is quenched, i.e., 
set irreversibly for a given realization of the polymer. 
We assume that there are no permanent dipole moments 
and all the dipoles are induced; thus is a quantity 

fluctuating about zero. 

For such a system, the Hamiltonian is given by 
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where the prime indicates that the terms with i = j 
and m = n are omitted from the first term. The above 
expression is an analogue of the polarizable field model 
proposed in Ref. [55H55] . developed here for the case of 


heteropolymers. The first term in Eq. (2.1) describes 
the electrostatic interaction between a pair of induced 
dipole moments, and the second term describes the self¬ 
polarization of the polymer, which can be regarded as 
the elastic energy of the local charge distortion [laiiTj. 
Here, G(ri,r 2 ) denotes the Green’s function, defined by 


G(ri,r2) 


1 

dTreeoki - ^21’ 
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where e is the relative permittivity of the dielectric 
medium the polymers inhabit. If we assume that wa¬ 
ter can be approximated as a uniform dielectric medium 
then e ~ 80 for the case of aqueous solutions at room 
temperature T = 293 K. 


A. Interaction between effective monopoles 


Here we make an interesting observation that the first 
term of Eq. (2.1) can be expressed in a form that de¬ 


scribes interactions between effective monopolar fluctu¬ 
ations, if (i) the polymers are semiflexible and (ii) the 
dipole fluctuations are parallel to the polymer back¬ 
bone. Taking the continuum limit, i.e., to —>■ s/a, where 
—1'/2 < s < is the continuous arc-length coordinate, 
these two conditions are summarized by the equation: 

p(s) = p(s)i(s) = p[s)R{s). (2.3) 

By defining the total polarization vector P{r), where 


Let us denote the position vector of the TO-th monomer 
(which can be, for instance, a nucleotide base-pair of a 


N 


P{r) = -J2 


(s)i^«(s)^(r-i^«(s)), (2.4) 
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We can rewrite the first term of Eq. (|2.1|) as follows: 


1 


V V im) ■ .„ 0 ) f 


mn 
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drdr' P[r)-Vr'G[r - r')-P{r'). (2.5) 


By partial integration, the expression (2.51 simplifies to 


drdr' Vr-P{r)Gir - r')Vr'■-P(r'), (2.6) 
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where V^. • P{r) is the polarization monopolar charge 
density: 


Wr-P{r) = - 




^ ^ r . 

= —^ / dsp^^\s)dsS{'i — R^^\s)) 


N r 


2 = 1 


= -I] P^^\--)S{r-R^^H--)) 

n ’ ’ 


-p«(-)d(r-i2W(-)) 

+ J dsp^^\s) 6 {r — R^^^s)) 


(2.7) 


There is a positive polarization charge of magnitude 
Pi{—^)la at i?i(—|), a negative polarization charge of 
Pi{—^)la at Ri{^), and a linear polarization charge den¬ 
sity of Pi{s)/a along the chain. Thus, for dipole fluctua¬ 
tions that are induced parallel to the polymer backbones, 
their interaction kernels can be transformed to monopo¬ 
lar Coulomb potentials, and their interaction can be re¬ 
garded as that between effective monopolar fluctuations. 


B. Quenched sequence-specific polarizabilities 

We now address the issue of sequence specificity for 
rigid polymers. In the continuum limit with m —>■ s/a 
(see above), a polymer i with a chosen sequence of 
monomers has also a quenched-in polarizability sequence 
a^*^(s), where a*'*^(s) of each monomer along the poly¬ 
mer chain i depends on its intrinsic physical properties. 
We need not specify at this point what the monomer po¬ 
larizability is due to or indeed how the monomers are 
defined. Assuming that the different monomers have a 
mean polarizability, which we call Q:q*\ then the differ¬ 
ences between monomers are reflected in the deviation 
(5 q:^*)(s) = a^''\s) — of the polarizability of each 
monomer segment from the mean polarizability value. 

In the absence of knowledge about how polarizabili¬ 
ties q;^*^(s) are distributed along each polymer i, we shall 
assume, for simplicity, that the values of the polarizabil¬ 
ities are drawn from a Gaussian probability distribution 


FIG. 1. Examples of sequences and sequence pairs that are 
studied in Sec. IV. Each colored circle represents a monomer 
segment, and distinct colors represent distinct polarizabilities. 
Sequences (a) and (b) are distinct, and their interaction free 
energy is studied in Sec. IV A, whereas sequences (a) and 
(c) are identical, and studied in Sec. IV B. Polarizabilities 
on sequences (a), (b) and (c) have zero intra-chain sequence 
correlation length, whilst those on sequence (d) have nonzero 
intra-chain sequence correlation length. The effect of nonzero 
correlation length is studied in Sec. IV C. 


! defined by its average and its variance. De¬ 
noting the sequence average over Pq, as ((...)), the average 
of polarizability is then 

((a«(s)))=ao, (2.8) 

and the sequence correlator is defined as 

(s) (s'))) = 9 ^a 6(3 - s')(5y. (2.9) 

On the other hand, if the monomeric polarizabilities 
along the polymer backbone are correlated over a nonzero 
length scale, we can approximate the sequence correlator 
by 

{{Sa^^\s)Sa^^\s')))c^g^afis-s')5,„ (2.10) 

where /(s — s') can be any decaying function of 
the argument such as a Gaussian, viz., /(s — s') = 

,_ ^ {a-s')2 

{\/^a)~^e 2 ^^, where tr is a typical length scale. 

Other choices are of course also possible. The form 
Eq. ( |2.9[ ) thus corresponds to a vanishing correlation 
length. 

In what follows we will be interested in the following 
three cases pertinent to two (i, j = 1, 2) interacting poly¬ 
mer chains [cf. Fig. [^: 

i) Polymers are distinct: q:^^^(s) =/= a^^^(s). There is 
a different probability distribution PQ,[ab)(s)] cor¬ 
responding to each polymer, and polarizabilities on 
distinct polymers are uncorrelated. Thus, we have 

{{a^^Hs)a^^Hs')))=al (2.11) 

ii) The pair of polymers are identical, i.e., a^^^(s) = 

(s) = q:(s). In this case, there is the same prob¬ 
ability distribution Pa[Q:(s)] of the polarizabilities 
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for either polymer. If we assume a vanishing corre¬ 
lation length, we obtain the following polarizability 
correlator on each polymer: 

((q;(s) a{s'))) = aQ+ g^aS{s — s'). ( 2 - 12 ) 


iii) The third case is similar to the second case (iden¬ 
tical polymers), except that the intra-chain corre¬ 
lation length of the sequence polarizability pattern 
is now nonzero. In this case, we have 

((a(s) a{s'))) =al+ g^af{s - s'), (2.13) 

where again /(s — s') can be any decaying function 
of the argument. 

Next, we will derive a field theoretical representation for 
the partition function of the system, which will enable 
us to study the behavior of the sequence-averaged vdW 
interaction free energy for the above three cases for a pair 
of rigid, rod-like polymers. 


III. FIELD THEORY FOR RIGID POLYMERS 
INTERACTING VIA INDUCED DIPOLES 


We are interested in evaluating the effective interaction 
(free) energy between a pair of rigid, polarizable poly¬ 
mers, where the dipole fluctuations have been integrated 
out. The task is made non-trivial by virtue of the fact 
that the interacting dipoles in the Hamiltonian (2.1) are 
coupled to one another. They can be decoupled by intro¬ 
ducing an auxiliary field ip and performing a Hubbard- 
Stratonovich transformation on the partition function. 
The resulting expression for the partition function is then 
amenable to a systematic approximation procedure using 
techniques of perturbation theory, which enables one to 
systematically approximate the effective polymer inter¬ 
action energy. 


A. Hubbard-Stratonovich transformation 


We can rewrite 


g 2 ) dR(')(m)dRO)(ri) ^ ' ' (3-1) 


_ g— ^ In det G 


Vipe 


8R(*)(m)8RO)(n) 


where 13 = and we recall that the prime denotes 

the exclusion of terms for which m = n on the same 
polymer. For a given chain conformation, but allowing 
for thermal fluctuations of dipole moments, the partition 


function is given by 

Z = n / dp<-'\m)e-^^ 

i,m 

_g-ilndetG J ^ J dr (V vf 
i,m 


xe ' 


(3.2) 

where we have made use of the Hubbard-Stratonovich 
transformation in Eq. ([33. Note that the dipoles are 
now decoupled, which enables a straightforward integra¬ 
tion over the dipole fluctuations. 

For polymers that are much more polarizable along 
their backbone than perpendicular to it, the induced 
dipoles prefer to fl uctuate parallel to the polymer back¬ 
bone [cf. Eq. (2.3)] in the direction of the local tangent 
vector of the polymer, t^*^(m). The dipole orientation 
is thus fixed by the polymer backbone orientation, and 
does not fluctuate if the polymer conformation is fixed. 
On the other hand, the magnitude of the dipole can still 
fluctuate in the whole accessible interval. Negative val¬ 
ues correspond to the dipole pointing in a direction oppo¬ 
site to the tangent vector, and positive values correspond 
to the dipole pointing in the same direction as the tan¬ 
gent vector. The functional integral over (m) then 
becomes a functional integral only over the magnitude 
(m). This involves evaluating the following product 
of Gaussian integrals: 


/ dp^'^{m) 




(3.3) 


i/S ^ (m )tW (m)- (m)) 


xe 


= Y 27rfcBTabl (to)- 


“f X] a^''^ {m).V(p(R^''\m))') 


The local polarizability of the polymer segments thus 
gives rise to a renormalization of the local dielectric con¬ 
stant. The partition function now reads 

Z = ZoJVif 

-i/3 

xe 3 (3.4) 

where we have defined the bare partition function Zq by 


Zq= Dip e 


— ^ f dr dr'ip(r)G ^ (r ,r')ip(r') _ — ^ In det G 


(3.5) 

and G ^(r,r') = — eeoV^5(i — r'), so that the Boltzmann 
average with respect to Zq is given by 

(. . .)o = Zg^JVip(. . _)e-§f‘irdr'<fi(r)G-'^(ry)v’(r')_ 


(3.6) 
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In Eq. (3.4) we have subtracted away the constants in¬ 


volving products of \/ ab)(77i). The meaning of the ex¬ 
ponent is as follows. Because there are no fixed charges, 
the mean electrostatic field is zero. On the other hand, 
dipole fluctuations give rise to electric field fluctuations 
iV(/7 in space. Owing to polarizability of the chains, the 
dielectric permittivity of a segment along a given chain 
is distinct from that of other segments, and also distinct 
from the dielectric permittivity of the surrounding space. 


Since the representation (3.4) entails a Gaussian inte¬ 


gral over auxiliary fields it can be evaluated explicitly, 
yielding 


In Z = In Z{ 
1 


0 


(3.7) 


--Indet ((eeoV^ -|-^cr^y(T')VaV6)5(r -r')) 


tions {i^W( m),given by 


F,eif[{H«(m),t«(m)}] (3.11a) 

i m^n 

d^GiR<^^'>{m), 9^G(fiW(TO), flW(^)) 

dRa^ {m)dRc^ (n) {m)dR^^^ (n) 

R^t[{R^^Hm),t<^^\m)}] (3.11b) 

= “ ^ (w)a(^) (n)tW (TO)tW {rn)t\?\n)t^l\n) 

mn 


i9^G(flW(TO),fi(^')(n)) d^G{R^^{m),R^^'i{n)) 

dRa^ {m)dRy^ (n) {m^dR^p (n) 


where we defined ^^{r) = / ds ta\s)tP{s)S{r — 

i?^®^(s)) in the continuum limit, where the subindices 
a,b = 1,2,3 denote the Cartesian components. The 
above expression for the partition function still retains 
the full dependance on the polarizability sequence (s). 


In deriving Eq. (3.10), there are singular terms that in¬ 


volve Coulomb Green’s function evaluated with two coin¬ 
cident field points on the same polymer, but their deriva¬ 
tives vanish and thus these terms do not contribute to Z. 


Using Eq. (2.2) and again going to the continuum rep¬ 


resentation, Tint can be re-written as 


B. Perturbation theory 

In order to proceed, we make use of perturbation the¬ 
ory. We expand Z to quadratic order in 

Z ~ ZpJVif e-ldfdreeo(Vvf 

x{l — i/3^a^*^(TO)(t^*^(TO) • V(p{R^^\m))p 

i,rrL 

+ X] X! a^^\n) 

ij mn 

X (t^*^(m) • V(p{R‘'^\m))p 

xp^\n) ■ Vv?(fi(^'Hn)))^} (3.8) 




k^T 


647r2(eeo)^ 




(3.12) 


i (*) (s) • t (J') (s') 3(i (') (s) • R,, )P^\s')- R,,) 


R^ 


R^ 


where Rij = ii*^*^(s) — i?^-l^(s'). The continuum limit of 
the interaction free energy is clearly non-pathological for 
different polymer chains (i ^ j). The above form of the 
interaction energy clearly reduces to the standard 
form for isotropic chains at large separation. 


Eollowing standard techniques of field theory [28] , we find 
the correlator 


(VW^(i^«(m))Vi^■V(i^^'Xn)))o 

, ^d^G{R<^^\m),R<-d\n)) 

= ’ —rr?-7-TT 

dRP{m)dRPin) 


(3.9) 


App lying Wick’s theorem to Z in Eq. (|3.8|) and using 

(3.10) 


Eq. (3.9), we obtain after re-exponentiation 

2’ ~ g —/5(T3olt-|-fint) 


Here, Fseii and Tint are respectively the free energy of 
individual polymers and the free energy of interaction 
between pairs of polymers (for given polymer conforma- 


IV. ROD-LIKE POLARIZABLE POLYMERS 


Let us specialize to the case of two rod-like polymers 
each of length i, and consider the three cases delin¬ 
eated above: (i) distinct polymers, (ii) identical polymers 
with zero intra-chain sequence polarizability correlation 
length, and (iii) identical polymers with non-zero intra¬ 
chain sequence polarizability correlation length. The 
P„[a«( s)]-averaged (or sequence-averaged) interaction 
energy for case (i) will be found to be identical to that ob¬ 
tained for polymers each having a uniform polarizability 
ao, whereas the sequence-averaged interaction energies 
for cases (ii) and (iii) will depend also on the variance 
around a^. 
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FIG. 2. A pair of rod-like polymers parallel to unit vectors 
and respectively and whose separation runs parallel 
to the unit vector r can equivalently be described in spherical 
coordinate geometry with r as the reference axis, 9i and 62 
as polar angles, and (j)i and (f >2 as azimuthal angles. The two 
polymers are co-planar if (f)i = (j> 2 - 


A. Distinct polymers 


After appropriate sequence averaging, the interaction 
free energy becomes 


{{Fint))=Fo (4.1) 

k^Ta^ 

3it(^\si)-Ri2Kt(^'>{s2)-Rl2) Y 

Fi 2 


Here, we have defined Fq as the sequence-averaged inter¬ 
action free energy of a pair of distinct rod-like polymers, 
a = ao/a is the average polarizability per unit length of 
each monomer, and we have let —^/2 < Si, S 2 < ^/2 p5] . 
The pair interaction free energy of distinct polymers is 
thus equivalent to that of polymers of uniform polariz¬ 
abilities Oo- _ 

Equation ( |4.1[) can also be expressed in the following 
form [see App. |A|: 


Fn = - 


ksTa 




fll 2 


dsi 


fij 2 


dS 2 


327r2(eeo)2 

6z{yir + {s2Z - Si)){y2r + (s2 - Siz)) 


\r + u\^ \r + u\° 

^9{yyr+J^z - Si))‘^iy2r + (s2 - Siz))^ 


\r u 
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(4.2) 


where r is the separation between the centers of mass 
of the two polymers, and u = S2t^‘^'^ — sii(i). The 
three parameters yi = • r, y 2 = • r, and z = 

f(i) . f( 2 ) characterize the orientational configuration of 
the pair of rod-like polymers, where r = r/r, and 
and t^ 2 ) are unit vectors characterizing the poly¬ 
mers’ orientations. The parameters j/i, j /2 and z are 
not completely independent, which we see as follows. 
Take r to be parallel to the reference axis in spheri¬ 
cal coordinate geometry, and write 6 for the polar an¬ 
gle and (j) for the azimuthal angle (see Fig. [^. We 


then have = (sin 0 i cos^i, sin 0 i sin(()i, cos 0 i), t(2) = 

(sin 02 cos (() 2 , sin 02 sin (^ 2 , cos 02 ), and f = (0,0,1). We 
obtain yi = ■ r = cos 0 i , y2 — ■ r = cos 02 , and 

Z = t(i) . t( 2 ) = where V" = 

cos ((/)2 ~ </>i)- The two rods are co-planar if 4 >i = 4 >‘ 2 -i- Qtt 
( or|V'| = l ;0 = O,±l,±2,...). 

In what follows, we study the behavior of Fq in two 
limiting regimes: (i) the near-held regime, r <g; where 
we consider two rods that are so near each other that 
the length of each rod can practically be assumed to be 
inhnite, and (ii) the far-held regime, r ^ 


1. Near-field regime 


We hrst consider the case where the length i of each 
rod is much greater than the separation between the rods, 
r i. The problem becomes analytically tractable if we 
consider the limiting case of inhnitely long rods. For 
such a case, the shortest length separation vector RI 2 
is perpendicular to both and As we show in 

App.lll Fq can be in that case approximated by 


M 2 a 2 kjiT z^ 1 

647r(eeo)2f2 7fr^|fi*2|4’ 


(4.3) 


where i ^^2 = r- + 

action energy is more attractive if the pair of polymers 
are more aligned (i.e., larger z), and vanishes if the two 
polymers are perpendicular to each other (i.e., z is zero). 
The result is invariant under z —)■ —z, which is to be 
expected as the polymers have inversion symmetry. 

The above form of the interaction energy is also fully 
consistent with the general form for the vdW interactions 
between two anisotropic rods in the non-retarded limit 
[30] if one takes into accou nt th e form of the anisotropy 
considered in deriving Fq. (4.3). 


2. Far-field regime 


Let us consider the far-held regime, i.e., the case of rods 
of length i which are short compared to their distance of 
separation, r ^ As the polymer lengths are hnite, we 
can rescale si, S 2 in units of £ so that — 1/2 < Si,S 2 < 
1/2. From Fq. (4.2), Fq can be put in a form where the 


integrand is dimensionless: 

^ M^alk^T FA 

° 327r2(eeo)^F V- 1/2 

2^ 62(^1 + (S2Z - Si)x){y2 + (52 - Siz)x) 


■ xuv 


jf-|- ccwj® 


^ 9 (yi + {S2Z - Si)xf{y2 + (S 2 - Siz)xY \ 


\r F XU 
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Here x = £/r, which is small in the far-held regime. Ex¬ 
panding in powers of x and integrating over si and S2 
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(i) 


{iii) 


\ 

\ 


(n) 


o 


(iv) 


least attractive when the orientations of the rods and the 
vector joining their centers of mass are all perpendicular 
to one another. For a pair of rods that lie in the same 
plane as the vector joining their centers of mass, as both 
rods rotate about their own backbones perpendicular to 
the plane through their centers of mass, the interaction 
free energy interpolates between power law decays of the 
form and r“®. 


B. Identical polymers 


FIG. 3 . Four configurations of a pair of polarizable rods: 
(i) \z\ = \yi\ = \y2\ = 1, (ii) \z\ = I,yi = y2 = 0, 
(iii) z = yi = 0 , |j/2| = 1 , and (iv) z = yi = y2 = Q. The 
interaction free energy is minimum for case (i) (where rods 
are perfectly aligned with each other and with the vector join¬ 
ing their centers of mass), higher for cases (ii) and (iii), and 
maximum for case (iv) (where the rods are perpendicular to 
each other, and each rod is also perpendicular with the vector 
joining their centers of mass). 


(see App. [§. we can put Eq. ( | 4 . 4 [ ) in the following form: 
r{( 2 : - 3yij/2)^ 


Fn = - 


M'^al k^T 
327r2(eeo)^r® 




yl - + yl)) 


+2/12/2(18- 35 (?/i +yl))z 




yl) - 2 )+ 


1 ^ 

^2 

-6 


}■ 


(4.5) 


The leading order term goes as which is character¬ 
istic of interactions between pairs of isotropic particles. 

Note that Fq is also invariant under the operations 
of (which changes the sign of both yi and 

z) and/or —>■ (which changes the sign of both 

?/2 and z). This is expected, as the free energy cannot 
depend on whether one defines one end or the other of a 
rod as its “head” (or “tail”). 

We next consider Fq for the following four illustrative 
cases (cf. Fig. [^: (i) |2;| = \yi\ = \y2\ = 1 , (ii) |z| = 
: 0, (iii) z 


1 , 2/1 = 2/2 = 

2/1 = 2/2 = 0 : 


= 2/1=0, IJ/2I = 1, and (iv) z = 


ksT /. I 4 t^ ^ 

327r2(eeo)^r® ^ “r ,^2 ^ 

M^gg fceT _ 3^\ 

32772(eeo)^r6 2r2 J 

kuT ( 3^2 

32772(eeo)+® V / 


(kl = I2/1I = I2/2I = 1 ), 

(|z| = l,yi = 2/2 = 0), 
(4.6) 

(j/i = z = 0 ,11/2I = 1 ), 


0 


(z = 2/1 = 2/2 = 0), 


to the order of 0 {r~^°). Fq is most negative for case (i), 
i.e., attraction between the rods is strongest if they are 
aligned with each other, and each rod is also aligned with 
the vector joining their centers of mass. The interaction 
free energy is least negative for case (iv), i.e., the rods are 


Next, we consider the situation where the polyme rs are 
identical, viz., a(^^(s) = a(^)(s) = a(s). Equation ( | 3 . 13 [ ) 
then becomes 


Eint 


ksT 

327r2(eeo)2 



ds' 

a 


a{s) a{s') 


[ Rh R!2 \ ’ ^ ^ ^ 

where we recall that R12 = r + The se¬ 

quence average is now given by 


((Ei„t)) = Eo + 5 F, (4.8) 


where ((...)) again denotes the average over the sequence 
probability distribution Pa and gives two contributions: 
the first, Eq, corresponds to the absence of correlations 
between the rods (which is the same quantity we calcu¬ 
lated for a pair of distinct polymers) 


En = - 




ds / ds' 


r3 

- 0-12 


327r2(eeo)2 

3(^(l)•i^l2)(^(2).i^l2)^" 


Ei 2 


( 4 . 9 ) 


while the second contribution 5 F corresponds to correla¬ 
tion of two identical polarizability sequences, assumed to 
be short ranged, of the form of a delta function [cf. also 
Eq. d^ ] 

■#( 1 ).#( 2 ) 

5 F = -, —” ' ■ ' * 


327r2(eeo)" 


ds 




3(i(i)-i2)(i(2).i?) 


( 4 . 10 ) 


where g = g/a is the root mean square value of 8 a = 
a(s) — pQ P6r unit length of each monomer, and Ft = 
_ ij(i)(s) = r + We have already 

calculated Eq in the foregoing subsection. We now cal¬ 
culate 5 F in the same two limiting regimes of near field 
and far field. 


1. Near-field regime 

Parallel to the near-field calculation for Fq in 
Sec. |IV A 1 [ the evaluation of the integrals in SF sim¬ 
plifies if we expand R about R* , the shortest length vec¬ 
tor that connects the two polymers at equal arc-lengths 


































FIG. 4 . Behavior of 7 = (vertical axis) for —1 < 

2/1) 2/2 < 1 (the two horizontal axes) for the following four 
values of 1/) = cos((()2 — 4 ‘i-) (the cosine of the difference in 
azimuthal angles of the polymers defined with respect to r as 
the reference axis): (i) 1)1 = — 1 , (ii) ip = 0 , (hi) ip = 0 . 5 , and 
(iv) Ip = 1 . We see that 0 < 7 < 1 , and 7 has a maximum 
value of 1 at 1/1 = —1/2 = ±1, and a minimum value of 0 at 
2/1 = 2/2 = ± 1 . 


(s = s'). Generally R* ^ Ri2'i equality obtains only if 
2/1 = 2/2 = 0 (in which case R* = R\2 = r). In the 
near-held regime, we can make the approximation that 
the rods are inhnitely long. Using Eq. ( 4 . , we hnd 

that 5 F is given by (see App. 0 


6 F = - 


3kBTMg'^xiyi,y2,z) 

16384-\/27r(eeo)^^r^ 


( 4 . 11 ) 


where the factor y(yi, 2 / 2 , 2 ) encodes the orientation de¬ 
pendence and is dehned by 

x(yi, 2/2, Z) = -^ + Uz + 41z2 

5(3-H0z-h3z^)(2/i-h2/2)^ 105(2/i + 2 / 2 )"^ ] 

1-7 4 ( 1 - 7 )^ j 

( — ' 1 ^ 

Here, 7 = 2 (i-z) ^ measure of the alignment or anti¬ 

alignment of a pair of identical rods. For —I < 2 / 1 , 2 / 2 ,2(’ ^ 
I [where we recall that ip = cos(( 2 ii 2 — (pi) is the cosine of 
the difference in azimuthal angles of the polymers dehned 
with respect to r as the reference axis], 7 varies between 
0 (for 2/1 = 2 / 2 ) and I (for 2/1 = — 2/2 = ±1 and z = — 1 ). 
This behavior is plotted in Fig. 

We make the following observations. First, we note 
that 6 F decays as r“®, which is distinct from the 
scaling of Fq in the near-held regime. The r~^ scaling is 


akin to a van der Waals interaction between an atom and 
a long thin rod. This is so because integrating along the 
lengths of two polymers subject to the constraint that 
s = s' is equivalent to an integration along the length of 
one polymer. 

Next, SF in Eq. ( 4 . 11 ) depends on z, in contrast to 
the z^-dependence of Fg. For polymers with uniform 
polarizabilities (which is what Fq effectively describes), 
the interaction energy does not differentiate between the 
“head” or “tail” of a polymer, i.e., there is inversion sym¬ 
metry (i.e., symmetry under the operation —)■ 

or <(2) ^ _t(2)^ 

or z —> —z). On the other hand, the 
breaking of inversion symmetry by SF rehects the seg¬ 
ment specihcity of the polymer polarizability. 

The behavior of SF is correspondingly enriched by the 
“head-tail” distinction of identical polymers. Compared 
to Fq for which a pair of parallel rods and a pair of anti¬ 
parallel rods are degenerate, the value of SF depends on 
7 whose value changes depending on whether rods are 
aligned or anti-aligned. Let us consider the behavior of 
SF for six orientation configurations. The first four are 
shown in Fig. whilst the last two are z = —l,2/i = 
1,2/2 —>■ —1 and z = —l,2/i =2/2=0 (which are the 
anti-parallel analogues of 2/1 = 2/2 = 1, 2 —1 and 2/1 = 
2/2 = 0 , z —>■ 1 , respectively). To study the behavior of 
SF for these configurations, we have to first determine 
the corresponding behavior of x{yi,y2, z) [ 5 T] : 


X(2/1,2/2,2)< 


164 

y/\ — Z 

= 336^2 
= 9 

, 36 

V2(1-7)6/2 

— 36 
“ 72 


(2/1 = 2/2 = 1,2 1) 

( 2/1 = y2 = 0,Z ^ 1) 

( 2/1 =2 = 0 , 12 / 2 ! = 1 ) 

(2 = 2/1 = 2/2 = 0) 

(2 = -1,2/1 = 1,2/2 -1) 

(2 = -1,2/1 =2/2 = 0) 


( 4 . 13 ) 


Applying Eq. ( 4 . 13 ) to Eq. ( 4 . 11 ), we see that whereas Fq 


was zero for perpendicular rods (z = 0 ), SF is negative 
and thus attractive for the same configuration. We also 
see that out of the first four configurations, SF is most at¬ 
tractive for perfect alignment (i.e., 2/1 = 2/2 = 2 = 1), and 
least attractive when the rods are completely orthogonal 
(i.e., 2/1 = 2/2 = 2 = 0 ). For sufficiently parallel polymers 
(z —>■ 1 ), X diverges as (1 — z)“^/^. Finally, each of the 
two anti-parallel configurations (z = —1,2/1 = 1,2/2 —>■ —1 
and z = —1,2/1 =2/2=0) is not degenerate with its cor¬ 
responding parallel counterpart (i.e., 2/1 = 2/2 = 1,2 —>■ 1 
or 2/1 =2/2 = 0 , z —)■ 1 ). Where x for the parallel-rod 
configurations diverges as (1 — z)“^/^ as z —)■ 1, the value 
of X for the analogous anti-parallel (i.e., z = —1) con¬ 
figurations can either diverge or not at all: in the case 
of z = —1,2/1 = 1,2/2 —>■ —1 (polymers sufficiently anti¬ 
parallel and the separation between their centers of mass 
is sufficiently parallel with the tangent vector of each 
polymer), x diverges, but with an intriguing (1 — 7)”®^^ 




















9 


0.0000 
- 0.0001 
- 0.0002 
-0.0003 
-0.0004 
-0.0005 

6 8 10 12 14 16 18 20 

(*) rja 

0 

-100 
-200 
-300 
-400 
-500 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

(n) '•/« 




FIG. 5. Comparison of the relative strengths of rescaled in¬ 
teraction free energies /o (solid line) and 5f (broken line) 
[defined in Eqs. (4.17l] for (i) g/ao = 3 and (ii) g/ao = 0.5, 
for each of the following four orientation conhgurations: (1) 
i /1 = 2 = 0, 1/2 = 1 and j/i = 1/2 = 2 = 0, (2) 1/1 = 1/2 = 2 = 1 
and 1/1 = 1/2 = 0, 2 = 1 , (3) i/i = 1/2 = 2 = 1 , (4) 1/1 = 1/2 = 0 
and 2 = 1 , (5) 1/1 = 2 = 0 and 1/2 = 1 and (6) t/i = 1/2 = 2 = 0. 
Here a is the monomer size. 


form as 7 —)■ 1 ; on the other hand, the sixth configuration 
(2 = — 1 , 1/1 = 1/2 = 0 ) has a finite and much lower value 
of 6 F compared to the second configuration. 

Thus at sufficiently short separations, “molecular 
recognition” of polymers will act in the following two 
senses: first, the r~^ scaling behavior of 6 F becomes suf¬ 
ficiently strong and dominates over the scaling of 
leading to pairs of identical polymers being selected 
over distinct ones; second, amongst the pairs of identical 
polymers that are selected, the ones that are sufficiently 
parallel {z « 1 ) or anti-parallel (2 « — 1 ) with 7 « 1 have 
the highest Boltzmann weight. 

What is the crossover length scale at which 6 F begins 
to dominate over Fg? Such a length scale can be found 
by comparing the relative strengths of 6 F and Fq. Using 
Eqs. (4.3) and (4.11), we have 


SF _ 3g'^T{yi,y2,z)£ 
Fo ~ 256V2Ma^r ’ 


(4.14) 


given by 


r(l/i,y2,2)= ^,^^'^5/,|9 + 142 + 4l2^ 

+ lOz + 3z^){yi + y 2 )^ 105 (i/i + 1/2)^ 1 

1-7 4(1-7)^ /■ 


Denoting by Rc the length scale at which SF/Fq ~ 1, we 
see from Eq. (4.14) that Rc is given by 


3g^'r(yi,i/2,2) ^ 

256V2Mal 


(4.16) 


The length scale Rc thus scales quadratically with g/ao- 
For rods that are perpendicular (z = 0), Rc diverges, 
and thus 6 F always dominates {{Fint})- The behavior is 
shown in Fig. where we define 


39^x{yi,y2,z)a^ 
^ ~ 256v^a§r5 


(4.17a) 

(4.17b) 


such that Fq = (ag/(647r/3(eeo)^-\/l — za^))fo and 5F = 
{al/ (647r/3(eeo)^\/l - 2 a^))5f. 

For orthogonal rods (2 = 0), Fq = 0 and thus ((Fint)) 
is always dominated by 5F. On the other hand, for rods 
that are sufficiently aligned 5F dominates over Fg at suf¬ 
ficiently short separations, with the crossover controlled 
by the ratio {g/afif'. To compare 5F for different orien¬ 
tation configurations, we should multiply the curves by 
a factor (1 — 2 )“^/^, from which we see that aligned rods 
are more attracted to one another than unaligned ones. 

To see what the crossover means in practice, let us con¬ 
sider a dsDNA rod that is ten base-pairs long (£ = 10a). 
For the 2 —1 configurations. Fig. then indicates that 
5F already starts to dominate over Fg at separations of 
order £ if the polarizability variance is three times the 
squared mean polarizability, but the crossover only oc¬ 
curs at separations of order Q.\£ (i.e., the monomer size) 
if the polarizability variance is half the squared mean po¬ 
larizability. Pairs of identical rods are thus more easily 
“selected” over distinct rods if the identical rods each 
have a more heterogeneous polarizability sequence. 


.8. Far-field regime 


Next, we consider the case where the separation dis¬ 
tance between the polymers is much larger than the 
length of each polymer, i.e., r ^ £. In this case, SF 
can be expanded in powers of £lr. Correct to quadratic 
order in £/r, we obtain (see App. H 


SF = — 


Mg^kBT 


h2 

12 


+ 0 {{£/r)f 


where the orientation-dependent factor T{yi,y 2 ,z) is 


327r2(eeo)2r® 


(4.18) 
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where the coefficients ho and h 2 are defined in Eqs. (E2). 


Erom the leading order term, we see that 5F is most 
attractive when the pair of identical polymers are aligned 
{z = — 2/1 = — 2/2 = 1 ) or anti-aligned (z = —l,yi = 
— 2/2 = 1 or z = — 1 , 2/1 = — 2/2 = —!)■ The leading order 
term is invariant under the inversion —>■ or 

f(2) ^ _t(2)^ 

or z —>■ —z, and the inversion symmetry is 
broken only at sub-leading order 0 ((f/r)^). 

To leading order, we then remain with 


kBTiM^al 

y-Cint// —- 


82 / 12 / 2 )^ 


+o{ie/rr 


327r2(eeo)^r® 

(4.19) 

i.e., the only effect of the identity of polymers is to renor¬ 
malize the overall sequence-averaged free energy. This 
seems again reasonable for the case considered, i.e., the 
far field regime. 


C. Identical polymers with non-zero sequence 
polarizability correlation length 


As in the previous case, we have = 

q;(s), while in addition the correlator decays over a 
nonzero length scale, Eq. (2.101. The interaction free en¬ 
ergy before sequence averaging is given by Eq. (4.7) and 
after sequence averaging we can express the free energy 
as the sum 


((Fint)) = Fo + SF, (4.20) 

where the average ((...)) again denotes the sequence av¬ 
erage, i.e., an average over the sequence probability dis¬ 
tribution Pa . Above, Fq is the part that does not involve 
sequence correlations, and 6F is the part that does in¬ 
volve such correlations. 

As before Fq has already been calculated [see case (i)]. 
The SF contribution is given in complete analogy with 
the case of short range correlations, except that 6{s—s') is 

,_ ^ (s-s')2 

substituted by a Gaussian f{s—s') = (v^cr)“^e 2,^2 . 

This leads to 

327r2(eeo)2 J_g/2 a 

W 2 “j 

In what follows, we explore the effects of having a 
nonzero sequence polarizability correlation length. We 
consider two situations: 1 ) large sequence correlation 
length (cr ^ i) and 2 ) short sequence correlation length 
(cr <C i). In the former situation, we find that Fq al¬ 
ways gives the dominant contribution to ((Tint))- In the 
latter situation, on the other hand, we will find the emer¬ 
gence of three qualitative scaling regimes for 6F: a far- 
field regime characterized by scaling, an intermediate 
near-field regime characterized by scaling, and an ex¬ 
treme near-field regime characterized by r~^ scaling. 


ds2 


1-1/2 a \p2/Ka 


(4.21) 


1. Large sequence correlation length, I a 


Let us first consider th e cas e of large sequence corre¬ 
lation length. Erom Eq. (4.21) we see that 5F becomes 
smaller as cr becomes larger, vanishing as tr —> 00 . The 
sequence-averaged interaction free energy of the identical 
rods then becomes no different from the interaction free 
energy between two uniformly polarized rods [i.e., Fq in 
Eq. (4.1)]. This is so, because a larger sequence corre¬ 
lation length brings about a stronger preference for the 
polarizabilities on each rod to be alike, and the rod thus 
appears to be more uniformly polarized. 

As before, we can consider the near-field {r i a) 
and far-field {i r a) behaviors. As i a, we 
can expand the Gaussian kernel ~ 1 — 

(si — S 2 )^/ 2 cr^. In the near-field regime, we consider for 
simplicity the case where 2/1 = 2/2 = 0 , and as before we 
approximate the rods by infinitely long ones. We find to 
subleading order 


SF = -= 


Mg^ksT 


1024-\/27r7r(eeo)^-\/l ~ + z)^/'^Lar‘^ 


X <{ 3(1 -b z -b 2z2 + 2z3) - z + z^)^\. (4.22) 


In the near-field regime, 5F has the same leading order 
r~^ scaling behavior as Fq, and the subleading co ntri- 
bution is of order r "^A comparison with Eq. (4.3) 


shows that the leading order contribution of 8F differs 
from Fq by a factor 


3(1 -b z -b 2z^ -b 2z^) 


1 - 


{yi-V2fi 

2{l-z) 




16-\/27rz^ 


Mapcr 


+0 


(4.23) 


Thus, in the near-field regime, the effect of dF is to renor¬ 
malize the prefactor of the vdW interaction energy, whilst 
the scaling behavior remains essentially unchanged. 
As Majl 00 , ((Eint)) becomes dominated by Fq. 

In the far-field regi me {£ r), the ratio x = ijr is 

small. Recalling Eq. (A2), R\2 = r + u where u = 
— and resc alin g si,S 2 ,cr — )■ sii, S 2 £, cri, we 

have [analogous to Eq. (4.4)] 


SF = - 


M^g^k^T 


r.1/2 
dsi 


1/2 

ds2- 




327r2(eeo)^ J- 1/2 J- 1/2 


6z(2/i -b {s2Z - si)x){y2 + (s2 - siz)a:) 


■ xu\' 


+ 


Jr-b xmJ® 

9(2/1 + {S2Z - Si)x)‘^{y2 + (S2 - Siz)xy 


XU 


10 


(4.24) 


Eor large correlation length, we can Taylor expand the 
Gaussian function to leading order. Then using Eq. (C2), 
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we expand 5F to quadratic order in x: 

Mg'^kBT 


- 


32-\/^7r2(eeo)^ or® 


1/2 ^1/2 

dsi / ds2 

-1/2 J- 1/2 


1 - 


(si - S 2 )' 

2o-2 


{fo + fiX + f2x'^)- (4-25) 


The coefRcients fi {i = 1,2,3) are defined in Eqs. (C3|. 


Evaluating the integrals over si and S 2 and restoring di¬ 
mensions to a yields 


SF = - 


Mg’^k^Ti 


32-\/^7r2(eeo)^ or® 


(4.26) 


^ { fo + 


Wi + W2 
12 


+0{MFa-'^r-^). 


e 

r2 


12 cr 2 


The quantities wi and W 2 are defined in Eqs. (C3). In the 
far-held regime, the contribution of 5F is always much 
smaller than that of Fq. The effect of the non-zero se¬ 
quence correlation length is to make 5F less negative. 


corresponds to a situation where the centers of the rods 
are directly opposite one another such that the vector 
joining these centers is perpendicular to the orientations 
of both rods; however, the rods are free to rotate in the 
same plane relative to each other. Dehning v = S 2 — si 
and u= {si + S2)/2, SF becomes 


dF = - 


ggksT 


e/2 .e ^ ^2 

du dve ^ 


32 v^ 7 r 2 (eeo)^CT J-e/2 J-e 

z 

(r2 -1- 2(1 — z)u2 -I- 4(1 -I- 

|(1 + z)2i;2-3(1-z)2w2 

(r2 -I- 2(1 — z)u'^ -I- |(1 -I- z)u2)®/2 


. (4.28) 


We can scale out the dependence on 1 — z by changing 
u —>■ uj— z) and v —>■ \f2vl\J\ F z\ 


SF = - 


gg kBT 


32-\/^7r2(eeo)^-\/l — z'^a 

p\f^t ..2 


(4.29) 


du 


dv e 


2. Short sequence correlation length, a <^ £ 

Next, we explore the consequences of having a short 
sequence correlation length. We can study three lim¬ 
iting cases: (i) intermediate near-held {a r <^ £), 
(ii) extreme near-held {r a i), and (hi) far-held 
(cr < ^ <C r). 

We begin with case (i), i.e., the intermediate near-field 
behavior. Let us rescale si, S 2 and o’ in units of r; equa¬ 


tion (4.211 becomes 


SF = - 


gg kBT 


fij2r 


fil2- 


\T — 


(■sj -s'lr 


dsi 


327 r 2 (eeo) 2 r 5 J_i/2r J-e/2' 
z 


dS2 


\/2/ko 


(1 -I- 2(522/2 - S 12 / 1 ) +sl + s^- 2zsiS2)®/2 
3 ( 2/1 -b 52Z - 5 i)( 2/2 -b 52 - SjZ) 

(1 -f 2 (s22/2 - S 12 / 1 ) + sf + si - 2 zsi52)®/2 j 

(4.27) 


.V 


As the dimensionless quantities Ijr 1 and a 1, 
we can approximate the upper and lower limits of the 
integral over u by 00 and — 00 , respectively, and approx- 

\-i / (S 1 -S 2 )'' 


imate the Gaussian kernel (v^^cr) " exp ( 


20-2 


by 


the Dirac delta-function (5(5i — 52 ). Physically this cor¬ 
responds to the fact that at separations large compared 
to the sequence correlation length, the monomer polariz¬ 
abilities on the same polymer appear to be uncorrelated. 
Performing the integration over 5i, we obtain the form of 
SF in Eq. (4.10), w hose ne ar-held behavior has already 
been studied in Sec. llVBTl and found to scale as r“'*. 

Next, let us turn to case (ii), i.e., the extreme near-field 
(r ^ a <^ £) behavior. Eor simplicity, let us consider 
the case 2/1 = 2/2 = 0 whilst keeping z arbitrary; this 


((1 -I- z)u^ — (1 — z)u‘^) 


(r^ -|- -|- 


(^2 


2)5/2 


After scaling u and v in units of a and dehning A = r/a, 
the equation above becomes 


SF = - 


gg kBT 


32 (ecQ) 2 


t 


du 


dve 1 +^ 


rr±ji (A^ 

„2 n 2 


K(l+^- (1 - Z)«2) 

(A^ -I- 


.w2)5/2 


.y2'^3/2 


(4.30) 


As cr/^ 1, we can replace the upper (lower) limit of 
each integral by 00 (— 00 ). For simplicity, we focus on 
the behavior for z —> 1. For small r/a, we hnd (see 
App .0 


SF -- 


249Mg2fcBT 


20487r®/2(eeo)^-\/l ~ z£ar'^ 

Used in conjunction with Eq. (|4.3[), we hnd 


(z^l). (4.31) 


SF 2A9V2g^e 


Fq 32-y^MQ;QtT 


(z^l). 


(4.32) 


Thus, SF has the same scaling behavior as Fq for r <C 
a I■ In contradistinction to SF as obtained for rods 


of zero sequence correlation length (see Sec. IVB 1), SF 
of rods of nonzero sequence correlation length exhibits 
a crossover from scaling to scaling at r < tr, 
and SF is larger than Fq by a factor proportional to 

{gyal){eiMa). 
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Finally, we come to case (ill), i.e., the far-field (cr <C 
£ <C r) behavior. To leading order in r, we can approxi¬ 
mate Eq. (|4.2ip by 


SF = - 


Mg^kBTcrjz - iyiy2f 
32-\/^7r2(eeo)^^r® 
.l/2a 

X / dsi 

J-l/2a J-i/2o 


(4.33) 


^/2o- , ,2 

ds2e 2 + 0{r ') 


-erf 


(£) 


2 e- 


Mg‘^kBTa{z - iyiy2y 
32-\/ 27r7r2 (eeo)^fr® 


- 2 


+ 0 (r-"). 


two sequences. We studied two situations: (i) the intra¬ 
chain correlation length a of the sequence distribution of 
polarizabilities is much larger than the polymer length £, 
and (ii) cr is much shorter than £. For the first situation, 
we found that the leading order contribution of 5F does 
not change the scaling behavior of the vdW interaction 
energy, and in the far-held regime the vdW interaction 
energy is always dominated by Fq. For the second situa¬ 
tion, three qualitatively distinct scaling regimes emerge. 
In the far-field regime (cr <C ^ «C r), we derive to the 
leading order 


6F = - 


Mg^kBT{z - 3^17/2)^ 
327r2(eeo)^r6 


(5.3) 


In the far-held regime, the only effect of having a nonzero 
sequence correlation length is a simple quantitative mod- 
ihcation to the prefactor of 5F, whose leading order 
scaling behavior and orientation dependence remain un¬ 
changed. At large separations, the sequence correlation 
length appears point-like if it is much shorter than the 
rod length. As expected, we recover the result Eq. (4.19) 
in the limit if a —>■ oo. 


V. SUMMARY AND DISCUSSION 


We have analyzed the van der Waals (vdW) interac¬ 
tion (free) energy between two rigid, rod-like heteropoly¬ 
mers with sequence specihc polarizabilities assuming that 
each monomer in the sequence has a different, hxed 
(quenched) polarizability. Two fundamentally disjoint 
cases were treated. In the first case, we assumed that the 
polymers are distinct, each with a different polarizability 
sequence distribution. In the second case, we assumed 
on the contrary that the interacting polymers are iden¬ 
tical, with the same polarizability sequence distribution. 
Surprisingly, the results for the interaction energy in the 
two cases are fundamentally different. 

As a summary of our results, we find that in the near- 
held regime (r <C F) for two distinct polymers, each hav¬ 
ing an average polarizability oq and M monomers, the 
interaction free energy has the form of Eq. (4.3): 


Fn = - 


/cbT z'^ 


1 


647 r(eeo) 2€2 a/I - z"^ |-Ri2l'‘ ’ 


(5.1) 


where RI 2 is the shortest length separation vector be¬ 
tween the rods, while in the far-held limit (r ^ £) the 
vdW interaction energy is in the form of Eq. (4.5): 


Fn = - 


M'^al kBT{z- ^yiy2f 
327r2(eeo)^r® 


(5.2) 


where r is the separation with yi = tFl-r/r, y 2 = 
and z = FF 

For identical polymers the interaction energy again de¬ 
couples into the Fq term, identical to above, and a correc¬ 
tion 6F due to the polarization correlation between the 


which obviously just renormalizes the magnitude of the 
interaction energy in the case of two distinct polymers. 
In the intermediate near-field regime (a r ^ €), SF 
exhibits r~^ scaling, with the form of Eq. (4.11): 


6F = - 


3Mg'^kBTx{.yi,y2,z) 

16384-\/27r(eeo)^^r^ 


(5.4) 


where xiyijy 2 , z) is given by Eq. (4.12): 


x{yi,y2,z) = 


1 


9-kUz + Alz^ 


Z (1 - 7 )®^^ 

5(3-b 13z)(yi-b 2 / 2 )^ , 105(?/i-b2/2)'‘ 


1-7 


4 ( 1 - 7 ) 


(5.5) 


and 7 = 2 (i-fzj • extreme near-field regime 

{r a i), SF exhibits r~^ scaling [see Eq. (4.31)], and 


is larger than fp by a factor proportional to g^ljMotf^a 
[see Eq. (4.32)]. Although the identity of sequences ap¬ 
pears to simply renormalize the strength of the vdW in¬ 
teraction in the regime of large separations, in the regime 
of small separations such an identity generates a distinct 
angular dependence in the dominant order, and the in¬ 
teraction free energy even acquires a novel scaling 
form in the interval a r <^ £. 

Without giving any calculational details, we also note 
that for a pair of Gaussian chains separated by a dis¬ 
tance larger than their radii of gyration and 
and which are characterizable by a Gaussian spatial dis¬ 
tribution of monomers and an isotropic orientation dis¬ 
tribution 


ti^^(s)ff^(s') = 0; (5.6) 

t^a^ (s) (s) = t^a'’ (s) ^ (s) = Sab (5.7) 

(where the overhead bars denote averaging over the 
monomer segment orientations in the high-temperature 
limit), the correlation of identical sequences only quan¬ 
titatively renormalizes the prefactor of the vdW in¬ 
teraction energy and does not introduce any essential 
change to its scaling behavior. At large separations, 
r [(7?g^^)^ -b , the vdW interaction energy 
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scales as whilst at shorter separations (of the order 
of the radii of gyration) it scales as r~^. 

We have demonstrated for a pair of rod-like polymers 
that the separation at which the correction due to 
polarization correlation, S-F, becomes comparable to Fq 
grows as (g/oo)^- At separations smaller than Re, SF in 
fact dominates the vdW interactions between them. One 
can interpret this as a sequence recognition effect, which 
would tend to separate interactions between pairs of iden¬ 
tical heteropolymers from interactions between pairs of 
distinct heteropolymers. We intend to explore the vari¬ 
ous subtle consequences and generalizations of this effect 
in the near future. Compared to the case of near-held 
rod-like heteropolymers the sequence recognition effects 
are much smaller in all the other explored situations. 


VI. ACKNOWLEDGMENTS 

BSL would like to thank G. Klimehitskaya and V. 
Mostepanenko for enlivening discussions. BSL and RP 
would like to acknowledge the hnancial support of the 
Agency for research and development of Slovenia under 
grant Nl-0019. 


Making use of Eq. ( |A2[ ) , we obtain 

R\2 R\2 

_ Z _ 3(yir -I- {S 2 Z - Si)){y 2 r F (s2 - si-z)) 

|r-I-|r-|-Mp 

(A9) 


Used in conjunction with Eq. (4.1), the above result leads 
to Eq. (lU^. 


Appendix B: Distinct rods — near-field calculation 

First, we show that the shortest length separation vec¬ 
tor i2*2 is perpendicular to both and . We de¬ 
termine the values and S 2 that specify the arc-length 
positions of the end-points of the vector R\ 2 , and cal¬ 
culate Eq using R \2 and shifted arc-length coordinates 
vi = Si — s* and z /2 = S 2 — S 2 . 

First let us minimize R 12 with respect to si and S 2 - 
From Eq. ( |A2[ ), we have 

-^12 =r^ + 2(522/2 - siyi)r -I- (sf -I- S 2 “ 2siS2z)- (Bl) 


Appendix A: Derivation of Fo, Eq. (4.1) 


From the minimization condition we obtain the following 
solutions s); and 


We denote the rod orientations by unit vectors and 

* ( 2/1 - y2z)r 

l-z2 ’ 

(B2) 

and the centers of masses of the rods by position 
vectors and 

* {yiz-y 2 )r 

'^2=1 2 ■ 

1 — 

(B3) 

R‘'^\si) = + sit^^\ (Ala) 

We then find that 


R‘'^\s2) = -1- 52t^^^. (Alb) 

Rl^=r + S2t^‘^'> - s*it<^^'> 

(B4) 


The vector linking one point si on the first rod to a point 
52 on the second rod is given by 


Ri 2 = R^‘^\s 2 ) - R^^\si) = r + u, 
where we have defined 


r = 
U = 


(A2) 


(A3) 

(A4) 


The orientation of the vector joining the two centers of 
mass is specified by r = r/r, and the orientational con¬ 
figuration of the pair of polymers can be specified by the 
three quantities: yi = 1)2 = and z = 

From these quantities, we deduce the following relations: 


fP) U = S2Z — 5 i, 


-U = S2 — 5iZ, 
u r = 52?/2 - siyi, 
u'^ = sl + s^ — 2 siS2Z. 


(A5) 

(A6) 

(A7) 

(AS) 


{yiz - y2)r ^(2) {yi - y2z)r ^(i) 

1 _ ^2 ^ 1 — z 2 ® 


and 


(W= 1 - 


yl+vl- 22 / 12 / 2 Z 

1 — 


(B5) 


It can be verified that R \2 ■ = .R 12 ’ = 0- Next, 

let us define shifted coordinates vi and V 2 '. 


Vi = Si - Si= Si 
V2 = S2 — S2 = 52 


( 2/1 - y2z)r 

1 -Z 2 ’ 

{yiz - 2 / 2 )?’ 

1 — z2 


(B6a) 

(B6b) 


We can express R 12 in terms of the shortest length vector 
and shifted coordinates: 


i ?12 — R\2 a ~ ■ 


(B7) 


expression for Fq: 


Using Eqs. (2.2), (4.1) and (B7|, we obtain the following 
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S27T^eeor J.ef2 J-e/2 { \Ri2\^ \Ri2\^ J 

k^Ta^ ^ 

327r2(eeor J-i/2-^-i^ -i/2-\ (|ifj2P + + i4^i2„y^zf 

Qz{zV 2 - l'l){v2 - ZVi) ^{ZV2- ViY{v 2- ZViY \ 

(|Jl*2p +^1 + ^2~ 2.ViV2Z)'^ (I-R 12 P + + ^^2 ~ 2viV2zY / 


(B 8 ) 


We now take £ ^ 00 . The second and third terms vanish 
upon integration, and only the first term yields a non¬ 
zero contribution. We thus obtain Eq. (4.3). 


Appendix C: Distinct rods — far-fleld calculation 

In the far-held regime, a; = £/r <C 1, and by rescaling 
s —> s£, we can make the following approximation. 


Ir-bitl r ^"|r-|-a:M| 


ux 


(Cl) 


~ r ^"(1 — 2 nr' 

+{2n{n + l){r-u)‘^ — 


For small x, the integrand {...} in Eq. (4.4) can be ex¬ 
panded: 


{■■■} — fo + fi X + f2 ■ (C2) 

Here, the coefficients of x" are given by /„ (n = 1,2, 3): 


/o = z^ - 6yiy2Z -b 9 y?y?, 

(C3a) 

fl = ViSi + V 2 S 2 , 

(C3b) 

f 2 = Wisf + W 2 sl + W 3 S 1 S 2 

(C3c) 


where 

vi = 6(3j/iy2 - z){{£,yl - l)y2 - 2yiz), (C3d) 
V 2 = 6 ( 3 j/i 2/2 - z)(2y2Z - {byl - l)yi) (C3e) 


TUi = 3[3(1 - 2byl + myi)yl 

- 42 /i 2 / 2 ( 35 y? - 9)z + 3(9?/? - l)z% (C3f) 
W2 = 3[3(1 - 25yl + 60y^)yl 

- 4 ?/i?/ 2 ( 35 ?/? - 9)z + 3(9?/? - l)z^], (C3g) 

W 3 = 3[4z^ - 76 ?/i?/22:^ 

-|-(2 — 22{yf + ?/|) -|- 310?/??/?)z 
-360y?yi + 12yiy2(5(y? + y?) - 1)]. (C3h) 


Integration over si and S 2 yields 


ri/2 

1 - 1/2 

. 1/2 

1 - 1/2 

.1/2 

1 - 1/2 


dsi 


dsi 


dsi 


. 1/2 

1 - 1/2 

. 1/2 

1 - 1/2 

. 1/2 

1 - 1/2 


ds 2 fo = {z -3yiy2Y, 


ds2 fi = 0 , 


dS2f2 


= 1^1 +y2- 50y?y? -f 60y?y|(y? -by?)] 
+yiy2[18- 35(y? -by?)]z 


+l[9{yl + yl)-2]z^. 


(C4) 

(C5) 

(C 6 ) 


Making use of the above results, Eq. (4.2) can be put in 
the form Eq. (4.5). 


Appendix D: Identical rods — near-fleld calculation 


Let us hrst derive the shortest length vector R*. The 
square separation is given by 

R^ = + 2{y2 - yi)rs+ 2{1 - z)s'^. (Dl) 


By minimizing R^ with respect to s, we obtain 


This gives 


yi - y2 

2(1-z) 


r. 


(D2) 


R* = r + 


(yi - y2)r 

2(1-z) 


(i(2)_t(l)), 


(D3) 


and 


(i?*)2 = (1 - 7 )r^ (D4) 

where 7 = ■ R* coincides with r (the vector join¬ 

ing the centers of mass of the polymers) if yi = y 2 = 0 . 
R* is not perpendicular to or but it is perpen- 
dicular to 

= i(yi+y 2 )r, (D5) 

_^*.(j(2) _^(i)) ^0. (D 6 ) 
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Let us express 


Appendix E: Identical rods — far-fleld calculation 




(D 7 ) 


where /x = s — s*. Because of the orthogonality of R* 
and we find 

= \R*\'^ + 2 {l- z)^l^ = (l- 7 )r 2 + 2 (l-z)/x 2 . (D8) 
Using the above results we have 


• R = ryi — (1 — z)s 
= ^(yi +J/2)r- (1 - z)n, 

tU) • R = r?/2 + (1 — z)s 

= ^{yi+y2)r + {i - z)fj., 


(D 9 a) 

(D 9 b) 


and 


fe/2 r^(i).^(2) 3(t(i).H)(t(2).ij)V 

L,f[-S> --1 


g /9_ (yi -1/2)^ 
2 ( 1-2 


d/i 


\ ((1 - 7 )r 2 + 2(1 - 2)^2)3/2 

|(yi + y2)V-3(i-z)Vl^ 

' ((l-7)r2 + 2(l-z)/x2)5/2 / ■ 


(DIO) 


Next we rescale y —>■ yj— z) and let £ —>■ oo, ob¬ 
taining 


rU 2 rt(i).i( 2 ) z{t^^').R)(t^'^).R) 

as ' 


l-t/2 


i ?3 




dfi 


2 > 2 


\/ 2 (i - 2 ) 7-00 I ((1 - 

i(yi +y2)^?’2 - f(i - z)ix 
((1 - 7 )r 2 -h/i2)5/2 
Stt 

512 y/ 2 (l — z){l — 7)5/2r5 

105(yi +^ 2 )^ ' 

4(1-7)2 


X < 9 + Ifc + 412= _ 5d±l“* + + »)’ 


1-7 


(Dll) 


Used together with Eq. ( | 4 . 10 [ ), we obtain 
SF = - 


iMg'^kBT 


16384 -^ 2(1 — z) 7 r(eeo)^(l — 

x2 

105(1/1 +y2)^' 


X 1 9 + 14x + 41x“ _ + »)’ 


1-7 


4(1 - 7)2 


(D 12 ) 


which is equivalent to Eq. ( 4 . 11 ). 


Using Eqs. (Dl) and (D 9 ), we can put Eq. ( 4 . 10 ) in 
the form 


dF = - 


knTgg 
327r2(eeo)^ 

X ds 


(El) 


-112 I ?’®(1 + 2 (y 2 - yi){slr) -f 2(1 - z){s/rYf 
Qz{yiy 2 r'^ + (1 - z){yi - y 2 )rs - (1 - zfs^) 

r 8 (l -h 2(1/2 - yi)(s/r) -b 2(1 - z){s/rYY 

9 ( 1 / 11 / 2 ^^ + (1 - z){]ji - y 2 )rs - (1 - zfs'^f \ 

+ rio(l + 2 ( 1/2 - yi)(s/r) + 2(1 - z)(s/r)2)io ]' 

Let us expand the integrand in the above expression in 
powers of small s/r: 


SF = - 


where 


kBTgg 


327r2(eeo)2r6 


ds 


^0 + ^1 (“) +^2 (-) |, 


/io = (3yiy2 - 2)^, (E 2 a) 

h\ = 6 {yi - y2){l + 10 yiy 2 - 5 z){ 3 yiy 2 - z), (E 2 b) 

/i2 = 3 ( 3 i/i - 12yii/2 + 120yii/2 -b 3 yj - SOOylyj 

+66O1/1I/2 + 120^1^2 - 1320 i/i 1/2 + 660 i/i 1/2 
+2z — SSyfz + 120 yiy 2 Z — 408yly2Z — iSy^z 
+876yly2Z — AOSyiy^z — 8z'^ + 63yfz^ 

— 164 i/ii/ 2 -z^ + 63 i/|z^-b 6z^). (E 2 c) 

Now, performing the integration over s, the term linear 
in s vanishes, and we obtain Eq. ( 4 . 18 ). 


Appendix F: Nonzero sequence correlation — 
extreme near-field calculation 


The evaluation of SF in Eq. ( 4 . 30 ) for z —1 is aided 


by the use of the following integrals: 


poo pOO 

/ du I dv 

/—00 


e i+- 


7 I-g 2 (l + z) 

-00" -b -b v'^Y 4(1 -b z )2 


(El) 


X <^ifn 


A2 


2(1 -b z) 


1 + z-A^ 


POO POO 

/ du I dv 

/ —00 


v'^e 1+*= 


A2 

57r3/2 


-K, 


A 2 


2(1 -b z) 


- t/( 5 -1 

(A 2 + 1 x 2 +-g 2)4 32 A 4 1 ^ 2 ’ ’1 + z 


l,4g 1 +^ 


poo poo 


1057r3/2 /5 A 2 

512A4 


where K,y{x) is a modified Bessel function of the second 
kind, and U{a,b,x) is a confluent hypergeometric func¬ 
tion. For small A we obtain the following leading order 
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contributions: 



e 


1,2 

1 +^ 


(A^ +v? + v"^Y 

^2 

(A^ +u^ + v'^y 

^2 


(A^ +v? + 


27 r(l + zY 
A4 


12A4’ 


TT 

32A4’ 


(F2) 
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